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RADIATION FROM A CURRENT FILAMENT 
DRIVEN BY A TRAVELING WAVE 


D. M. Le Vine and Robert Meneghini 

Goddard Space Flight Center 
Greenbelt, Maryland 


INTRODUCTION 

This document presents a solution for the electric field radiated from an arbitrarily oriented 
current filament located above a perfectly conducting plane and driven by a traveling current 
wave. Both an exact and an approximate solution are obtained. The approximate solution 
is obtained by simplifying the integrand in expressions for the magnetic vector potential 
in a conventional manner (i.e., the “far-field” approximation) and, as a result, yields only 
the radiation terms of the electric field measured by an observer in the fraunhofer region of 
the filament. On the other hand, a closed-form expression for the fields can be obtained 
without resorting to approximation in the special case in which the current wave propagates 
along the filament at the speed of light in the medium surrounding the filament. This “exact” 
solution predicts both radiation and near (e.g., induction) fields and applies regardless of 
the distance of the observer from the filament, length of the filament, or frequency. Both 
solutions apply regardless of the current waveform: the only restriction on the current wave- 
form required in the approximate solution is that it be a propagating wave, and the only addi- 
tional restriction on the current waveform required for the “exact” solution is that the wave 
propagate with the speed of light in the medium in which the filament resides. 

This problem arose during a study of radiation from lightning. The signal radiated by a light- 
ning flash (called an “atmospheric” or just “sferic”) is important as a noise source in the 
design of communications systems and is also as a diagnostic in the study of lightning itself. 

It may also be possible to relate the radio-frequency radiations from lightning to meteorolog- 
ical observables of the parent storm cell: for example, it has been suggested (References 1 , 

2, and 3) that a high frequency of occurrence of lightning is related to the existence of 
tornadoes. 

The physical nature of lightning permits analysis at radio frequencies, in terms of a piecewise 
linear model of radiating electrical elements. For example, a typical cloud-to-ground lightning 
event begins as a series of ionizing steps (stepped leader) along which current propagates to 
establish a conducting channel between the charge in the cloud and the ground. The leader 
process is followed by a large current surge (the return stroke) which exchanges charge between 
the cloud and the ground. The current in the return stroke propagates along the channel at 


velocities approaching the speed of hght (Reference 4), and the shape of the current wave 
is roughly exponential (Appendix A). 

An obvious natural beginning point for the calculation of the fields radiated by a piecewise 
linear model is to use the ‘Tar-field” approximation to calculate the fields radiated by each 
current element in the model. Unfortunately, questions of validity arise when this approxi- 
mation is used for lightning, because interest often Ues in relatively low-frequency behavior 
(the peak energy radiated by a return stroke is centered near 10 kHz) and because the light- 
ning channel is often long enough (on the order of a few kilometers) so that being in the 
far field cannot be guaranteed even at relatively large distances from the flash. On the other 
hand, an exact solution exists for this problem whenever the propagating current wave 
travels at the speed of light in the medium surrounding the current element. Although not 
precisely appropriate to hghtning because of the restriction on the velocity of propagation, 
this solution is close enough to reality to provide insight into the near-field and low-frequency 
behavior of the field electromagnetic fields; it also provides a standard with which to assess 
the applicability of the far-field approximation. 

The purpose of this document is to present both the exact and the far-field solutions for 
representing radiation from lightning and to use these solutions to illustrate the effect of 
various channel parameters on the spectra (i.e., magnitude of the Fourier transform) of the 
fields radiated from lightning. The following section presents the exact solution for both 
the electric and magnetic fields radiated from a current filament driven by a traveling wave 
of arbitrary shape. The next section presents the equivalent solution, obtained by making 
far-field approximations, and the following section compares the two solutions. This com- 
parison indicates that the approximate solution fails at low frequencies in a regime which 
may be of interest in lightning studies. The comparison is followed by examples of the effect 
of filament orientation, velocity of propagation, filament length, and the observer’s orienta- 
tion on the spectra observed from a typical return stroke. Finally, a model is presented for 
a transfer function for the radiation process. 

EXACT SOLUTION 

Consider a current filament, J(r, t), above a perfectly conducting plane (the z = 0 plane) as 
shown in figure 1 . Assume that J(f, t) has the form: 

j(7,t) = ■^f ^t- (1) 

where £ is a unit vector parallel to the filament, and in the direction of current flow,T is the 
position vector of a point on the filament, and c = 3 x 10^ speed of light. Thus, 

J(F,t) is a current pulse with shape fit) which propagates in the £ direction at the speed of light. 

A solution for the Fourier transform,* of the radiated electric field, E(r,t), is to be 

obtained here by Fourier transforming Maxwell’s equations (with respect to time) and then 
solving the transformed equations. This, in turn, is done most conveniently in terms of a 


*The superscript “tilda” denotes a Fourier transform. For example: g'(i^) = y* g (t) exp (j2nvt) dt. 
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magnetic vector potential, A(r, i^). By noting that J(T, »») = fi f (i^)expO'k $ ' r) where k = 2 w/c 
and employing the image theorem for currents, the following result is obtained for the i* 
(rectangular) component of the magnetic vector potential: 


Aj (t,v) = fi £. J' f (v) ej^ ® * *■ 


JkR 


47tR 


— e, 


jjkRj 


‘ 4jtR, 


dr' 


where the integration is over the filament and where 
R 


R. = 


^0 “ ~ (yo “ y')^ + ~ 

(x - x')^ + (y - y')^ + (z + z')^ 


'^0 ~ ^1 


( 2 ) 


(3a) 

(3b) 


e. = 


+1 if i = X or y 
-1 if i = z 


(3c) 

(3d) 


rg = the position vector of the observer 

By transforming to a coordinate system in which one axis is parallel to the filament, equation 
2 can be reduced to a single integral. Let this axis be the z-axis, and denote coordinates of 
the observer and filament expressed in the new system byf' and r", respectively: 


Aj (r g , v) = M Cj 7 (v) 


/ 


Jkz' 


ejkr 

4irr 


— e. 


e^ 

4nr, 


dz' 


(4) 


where 


r =-^(Xg - x”)2 + (y'g - y")2 + (z^ _ z'')2 

^■i = + 4 [ Zg cos 0 - Xg sin 1 [ z" cos <t> — x'' sin <}> ] 


(5a) 

(5b) 
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and where x" and y'' are the x- and y-coordinates, respectively, of the filament as seen in the 
rotated reference frame. The rotation was accomplished by first rotating counterclockwise 
about the z-axis through the angle 6 until the x-axis lay in the plane formed by C and the 
z-axis, and then rotating counterclockwise about the new y-axis through the angle <p until 
the z-axis was parallel to C. (See Appendix B.) 

Now, in the first integral in equation 4 (i.e., the integral over the source terms), make the 
change of variables: 


U = r + (z” — 2 q) (6a) 

and in the second integral let 

Uj = rj + [(z" — Zq) + 2 (Zq cos^ 0 — Xq sin 0 cos 0)] (6b) 

This results in the following: 


Aj (r Q , v) = M Cj f (p) 




47TU 


du — 


(7) 


— €■ exp [jk (— Zq +2 sin 0 (Zq sin 0 + Xq cos 0 ))] 



Uj(a) 


ejkui 

4;rUj 



By transforming back to the original reference frame (Appendix B), and noting that 



A, A A A A 

where £ = C - 2(£ • z) z is a unit vector in the direction opposite to the image current, the 
following vector form is obtained for the magnetic vector potential: 
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■% 


A ( r. , I') = ixf(v) 


■Uflj) Jku 


• Fg y" 

( U(a) 

- i' • "o /" ** ^ du'l 

U'(a) 47TU ) 


du 


( 8 ) 


where for i = a, b; 


and 


U(i) = IPi 1 

- t • Pj 

(9a) 

u'(i) = Ip! 

1 - • p[ 

(9b) 


Pi = h 
Pi = ^0 


— r. 


r[ 


(9c) 

(9d) 


where r and r. are the position vectors of the end points of the current filament, r is the 
position vector of the observer, and = r^ - 2(r^ • z)z and r[^ = r^^ - 2(r,^ • z)z are the positi 


a ^ a 

vectors of the end points of the image current filament. 


position 


The electric-field intensity, E(r^ , 1 ^), can be computed from equation 8 by means of the 
formula: 


E ( , I.) = jkc 


A(f,i^)+ — V (V-A) 
° k^ 


( 10 ) 


The procedure is relatively straightforward and can eventually be written in the following 
form: 


E(r^,i^) = yjme7(y) 


4 7rp/ 


4 7Tp. 


e(Ph) - 


4 7T p^ ^ ^ 


, , eikpl 


( 11 ) 
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where for i = a,b: 


e(p.) = 


1 + g • V Pj 
1 - ^ • V Pi 


e (Pi) + j 


kPi 


e'Cpf) = 


1 + ^' • vpj 

1 - • VPj' 


^(Pi) 


^'(Pi') = 


^ - (^ • VPi) VPi 
Vl - (t • VPi)^ 

- (t' • VPpVPi 


A 


e'(p[) + j 


VPi 

kp; 


i = Vk - 


= ifo -fi' 


<yo - = 


■fo I 


(12a) 

(12b) 

(12c) 

(12d) 

(I2e) 

(12f) 


Note that vectors^ (p.) and ^ (p.') are unit vectors in the direction of the radiation component 
of electric field associated with the source and image current, respectively. Note, also, that 
the solution contains four terms, each of which involves the coordinates of the observer and 
one of the end points of the source or image filament. The two terms associated with e (p.) 
involve the end points of the source, and the two terms associated with ^ (p.') involve the 
end points of the image. Consequently, the radiation from the filament can be treated as if 
it were emanating from the end points of the source and image filament only. 

In the solution for E(F^ ,v), note the presence of the terms that are dependent on Vp./kp.. 
These terms represent inductive (i.e., nonradiation) components of the electric-field intensity 
and become negligible in the limit of large kp.. 

The magnetic-field intensity associated with the current filament can also be obtained from 
equation 10 by means of the formula: 

H(r„,»;) =— V X f (?„,»') (13) 

Op" 
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The calculation is again straightforward, although somewhat lengthy, and may be siroplified 
with the aid of the vector identities: 


VX (^0) = V0 X ^ X ^ 


(14a) 


v(0 0) = <j)Wd + e v<i> 


A 


(14b) 


Note that, for i = a,b, VU(i) = Vp. - £ and similarly that, in the case of the image, VU (i) = 
Vpf -^;thus: 


H(r„,i;) = f(j^) 


jjhPjj _ — 

h(Pb) - — - h(pj 


4irp^ 


4 7TP, 


(15) 


7 h'(p;) - — , h'(p:) 


4^Pb 


4 7tp' 


where for i = a,b: 


A 

h(Pi) 


h'(p;) 


CD 

»> 

1 +'2 ’ VPi 

V 1 - t • VPj 

h(Pi) 

gjk ^ • r j 

I 1 + ^'* Vp| 

\ 1 - • Vpj 

^ ' (P|) 

VP: X ^ 




V 1 - (^ • VPj)2 
vp[ X 

Vi - & ' • VPi')^ 


(16a) 


(16b) 


(16c) 


(16d) 


A. 


The unit vectors h(p. ) and h'(p.' ) are in the direction of the magnetic-field intensity associated 
with the source and image current filament, respectively. 

Note that the solution for the magnetic intensity does not contain terms which depend on 
powers of 1/kp. higher than the first. Also, the electric and magnetic fields associated with 
each end point (i.e., those attributable to one of the four terms in the solution) are orthogonal. 
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and the fields associated with the source filament only or the image filament only are orthog- 
onal, but the total electric and magnetic fields are orthogonal only in the special case in 
which the image and source are coplanar. That is, in the absence of boundaries, the total 
radiation fields from a single source are orthogonal, but the boundary, acting as an effective 
second source, can result in nonorthogonal electric and magnetic fields, even in the far field. 

Finally, the magnitude of the radiation components of the electric and magnetic fields asso- 
ciated with a given end point are proportional, the proportionality constant being the char- 
acteristic impedance of the medium, . As previously stated, these fields are also orthog- 

onal. Therefore, in some respects, the fields associated with each end point behave as if 
there were a small current source located at the end point. 

APPROXIMATE SOLUTION 

A solution is presented here for the problem treated above, radiation from a current filament 
driven by a traveling wave and located above a conducting plane, by employing the far-field 
approximations; that is, assuming that the observer is far enough from the current filament 
that kp » 1 and kL^ /p « 1 , where L is the length of the filament and p is the distance 
from the filament to the observer. When this approximation is made, the restriction to 
currents which travel at the speed of light, c, is longer necessary and the following current 
waveform which propagates with speed v in the C direction can be treated: 


X • T" 

J(7,t) = ^f(t —) (17) 

V 

Other than this minor change, the solution is formulated initially as in the proceeding section. 
Thus, the magnetic vector potential is: 


A(r^,t') = 


f 


filament 


JkR 


4 7tR 


A 

£ - 


4 7tR 




dr 


(18) 


where all quantities are as defined in equation 3. 

The electric field is^btained from the preceeding equation by means of equation 1 0, and 
noting that V[ V* (£0)] = (£• V) V 0 for the constant vector £, the following result is obtained: 
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E ( , i^) = jk yfiiie f (v) J e- 

filament 


/ e»,.. [^] 


$ + — ($ • V )V 
k2 


ejkR 

47tR 


^ ' + — (1^' • V ) V 
k2 


4n R, 


( 19 ) 


To evaluate the integrals in equation 19, assume that the distance from the filament to the 
observer is larger than the length of the filament and expand R and Rj in a power series about 
the center of the filament. Thus, in the case of the source filament. 


R = 1 r - r ' I = I ( r„ - ) + ( r^ — r ') I 

= V'^0 - Tc + K - r ' + 2(r^ - r^) • ( r ^ - f') 

1 Pre-F'!^' 

I r — I + I — r ' 1 cos 0 + — — 

^ ^ " 2 1^ _ . I 


( 20 ) 


where 


cos 0 = 


( •'o - ^c) • ( re - *■ '> 


If - r I It - r ' I 

o c c 


( 21 ) 


and where r^ is the position vector of the filament center. Now, let = Ir^ - r^ I, and assume 
that 


kL^/p^ « 1 (22a) 

kL » 1 (22b) 

where L is the length of the filament. Then, R = p^ + I7^ -T' I cos 0 for use in the exponential, 
and R = p^ elsewhere. Repeating this argument for the image distance Rj and substituting 
the results into equation 19, it follows that: 
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!(»')- 


B(t^,u) = jk-^tije 7(v) 



ei^Pc 

4irp^ 


JL 





4 7tp; 



( 23 ) 


where: 




Ir^ - r ' 1 

I(^) = J Qi2lTU |_ y J 

ej27rp 

c cos (j> 

filament 


_ _ 


df ' 


(24a) 


I' (i') 


/ 

filament 



ej27ri^ 



dr ' (24b) 


The vectors associated with the image filament are; 


- 2 (T, 

A . A 

• z ) z 

(25a) 

f| = r' - 2(f' 

A . A 

• z ) z 

(25b) 

P'r = Ir. - 1 

0 c 


(25c) 


To evaluate integrals Ifi") and r(t»), rotate to a coordinate system in which the filament is 
parallel to the z-axis as in the preceeding section and as described in Appendix B , and then 
perform the z-integration in a coordinate system with z-origin at the filament center. Note 
that, in this coordinate system, Ir^ -T' I cos (p = (z^ - z") £• Vp^ ; therefore 
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Uy ) = f [t] e~j 2 Wc[(z" - z^) ^ • Vp^] dz” 


= f ’ VpJ z" d2" 

-L/2 


= I L sine [ */^ kL (tj — ^ • V p^) ] | 


where 77 = c/i^. Similarly, in the case of the image, 


I' (v) = I L sine [ 54 kL( 7 ? — £' • Vp' )]| 


On transforming back to the original reference frame, replace z and z' by z, = C • r„ and 

A C C C C 

z' =£' - 7 . 

C c 

The remaining vector operations are straightforward and yield the following results: 


r t + — (^ • V ) V 

L 1,2 


1— - hr 

J I [ 


(£ • Vp^lVPj, 


1 + j 


(£ • Vp )Vp^ 




-L 


'1 


1 + j 


• vp;)vp; 
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Now, keeping only terms of lowest order in kp^ , the following results for electric field are 
obtained: 


! jk 


I(»^) - 


(28) 


- (^ • vp;)vp;j I 


where 


— 


I(j') = jkLe+jl^^ * rc)sinc [ '/2 kL (77 - t • Vp^) ] 


(29a) 


I’ (^) - jkLe"''J^ ^ sine [ V2 kL (tj — 


vp;)] 


(29b) 


The magnetic field intensity can be obtained from equation 1 3 in a straightforward manner. 
After making the far-field approximations and keeping only lowest order terms in kp^, the 
following form results: 


H(r;,i^) = f(v) 


(Vp^ X h 


eikPe 

4 ttPc 


I(t^) - (Vp; X t') 


4 7tp; 



(30) 


The solution in equation 30 is valid for the arbitrary index of refraction, 7 ?; however, in the 
special case, 77 = 1 , which is the assumption used above in the treatment of the “Exact 
Solution,” the far-field solution can be put into a form even more similar to that obtained 
for the complete solution. Adopting the notation used in the preceeding section, the follow- 
ing forms are obtained for the far-field solutions when 77 = 1 : 


E ( r^ , j') = 2j V P/^ f (*") 


4 7tp^ 




- ^ • Vp^)l - 


A ^ R ^ ^ 

4 7tp^ ^ <= 



(31a) 
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VPo) 1 - 


H(r„,i^) = liliv) 


eJkPc 

4 7TP^ 


h (p^) sin 


i/zkLd - $ • 


eJ’^^c - 

^ h'(p')sin [y 2 kL(l 

4 7TP^ 





(31b) 


where all quantities are as previously defined except that (p^) and e (p' ) are the vector 
components of equations 12a and 12b associated with ^(p) and (p'), respectively; that is, 
associated with the radiation terms only. 

Equations 31a and 31b obviously have a form similar to the radiation component of the 
complete solution, except that, whereas the complete solution (equations 1 1 and 15) gave 
the appearance that the radiation came from the ends of the filament, the far-field solution 
appears to emanate from the filament center and, in comparison with the complete solution, 
is weighted by a spatial ‘‘transfer function,” 2j sin [!/ 2 kL(l - £ • Vp)] , associated with the 
filament. 


COMPARISON OF EXACT AND APPROXIMATE SOLUTIONS 

The solution presented in “Exact Solution” is exact in the sense that no mathematical 
approximations were needed to obtain a solution in closed form. For this reason, this solu- 
tion provides a useful standard with which to compare the approximate solution obtained 
above. Comparison of the two solutions yields insight into the usefulness of the far-field 
approximation and the nature of its shortcomings. 

As a first point, note that making the far-field approximation in the exact solution (equations 
1 1 and 15) yields a result identical to the approximate (far-field approximation) solution 
obtained above (equations 31a and 31b). The details are given in Appendix C. In making 
the far-field approximation in the exact solution, the approximations must be made after 
the integrations over the filament have been performed, whereas, in deriving the approx- 
imate solutions (as in “Approximate Solution”), the same approximations are made before 
performing the integrations. Appendix C shows that the end result is the same, giving 
credence to the approximate solution. 

The conditions imposed by the far-field approximations used to obtain the approximate 
solution are that, in addition to L/p « 1 , kp » 1 and kL^/p « 1 which would suggest, 
for fixed p and L, a frequency band in which the approximate solution faithfully represents 
the true solution. Certainly, as kp becomes smaller, the nonradiation terms in the exact 
solution become significant, and the approximate solution, which contains only radiation 
fields, must fail. Appendix D contains a more quantitative assessment for an observer 
located on the surface. It shows that, at low frequencies, as expected, the far-field solution 
fails to agree with the exact solution whenever the nonradiation terms become significant. 
The exact solution is actually singular in kp, but the far-field solution is not. However, 
the solutions do not agree at low frequencies even when only radiation terms are compared. 
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Appendix D shows that the power spectrum for the radiation terms of the electric field as 
predicted by the exact solution becomes a constant times the spectrum of the current wave- 
form, I f(i^) I as k -»■ 0, whereas the far-field approximation is proportional to k times | fii') | 
as k 0, so that even the radiation terms exhibit a different low-frequency behavior. 

The relative low-frequency behavior of the exact and far-field solutions is illustrated in 
figures 2 through 6 for a current waveform with exponential shape. (A three-exponential 
model is used for the current. This waveform is shown in figures 7 and 8 and is discussed ' 
in Appendix A. The spectrum, I f(»^) I, associated with this waveform is constant at low 
frequencies and decreases as Ijv for large frequencies.) In each of figures 2 throu^ 6, 20 
times the logarithm (base 1 0) of the magnitude of the electric field is plotted along the 
ordinate, and frequency is plotted along the abscissa. In figures 2 through 4, the observer 
is located on the surface at increasing distance (5, 50, and 500 km, respectively) from the 
origin. The current element is a vertical filament 1.5 km long with one end at the origin 
(i.e., on the surface). Each figure shows three curves, which represent the exact solution, 
the far-field (approximate) solution, and the radiation component only of the exact solution. 
Because of the nonradiation terms, the exact solution appears as iTfi') l/k near k = 0. Because 
I T(p) I is a constant for k = 0 for the exponential waveform, this solution is singular at zero 
frequency. This singularity is evident in each figure where, for low enough frequency, the 
exact solution begins to increase with decreasing frequency. The asymptote for this low- 
frequency behavior is a straight line with (negative) slope of 20 dB/decade. On the other 
hand, the radiation component of the exact solution approaches a constant times 1 7(1^)! for 
small as is evident in each of the three figures. The low-frequency behavior of the approx- 
imate solution is different from both of these asymptotes; it decreases as 1/u for low fre- 
quency. The effect of increasing the distance from the observer to the source is both to 
decrease the magnitude of the spectrum at a given frequency and to shift toward lower fre- 
quency the point at which the three curves— exact, radiation, and approximate solutions— 
begin to differ. This is because the criterion kp » 1 is dominant for this case. The effect 
of decreasing the element length is illustrated in figures 4 through 6 in which the element 
length is decreased from 1 500 m to 1 50 m to 15 m, respectively. 

In the limit of very high frequency (high enough so that only the radiation terms in the exact 
solution need be considered), both solutions predict power spectra which are proportional 
to I f(p) I times an oscillating factor. In the exact solution, this factor oscillates between 
finite upper and lower limits determined by the geometry; in the far-field solution, this 
oscillation is sin [VakLfl - 2 • Vp^ )] and, in particular, has nulls that are not usually present 
in the exact solution. On the other hand, the peak values predicted by the two solutions are 
essentially the same whenever the length of the filament is short enough so that p^ — Pb • 
Therefore, if only the upper envelope of the magnitude of the power spectrum is required, 
it is proportional in both cases to I f(t») I, and the criterion that the proportionality constants 
be equal is that, for i = a,b: 
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Pi 


^(p^) • ^ (32) 



which is independent of frequency. This requirement will be satisfied if L/p^, « 1 . There- 
fore, although the requirement kL^ /p « 1 would indicate a failure of the far-field solution 
at large k, this does not happen if only the upper envelope of the power spectrum is required. 
In this case, the purely geometric criterion given in equation 32 is appropriate. 

The high-frequency behavior of the exact and approximate solutions are illustrated in figures 
9 and 10. Again the three-exponential waveform (Appendix A) has been used, and the ob- 
server is located on the surface at 10 km from the origin. In figure 9, the filament is vertical, 

1 .5 km long with one end at the origin. In figure 10, the filament is horizontal with its 
center 4 km above the surface at the origin and oriented so that the filament and observer 
are coplanar. Note that the exact and approximate solutions are essentially identical for 
frequencies above about 1 kHz (kp = 1 ). (Nulls in the approximate solution do not show in 
the plot because of computer sampling and roundoff errors.) The different slope of the spec- 
tra for the horizontal and vertical elements is attributable to the factor, L(1 - £ • Vp)— the 
effective length of the filament which is essentially L for the vertical element, but is consid- 
erably smaller for the horizontal element. Note that, for the vertical element, no apparent 
difference exists between the exact and approximate solutions (except near the nulls), even at 
frequencies above 1 MHz where the requirement, kL^ /p « 1 , is no longer valid. 

EXAMPLES OF SPECTRA 

The actual spectrum measured by the observer depends on the orientation of the current 
element, the polarization of the observer’s receiver, and the height of the observer above 
the surface, as well as the physical properties of the filament (i.e., current waveform, velocity 
of propagation, and element length). This section contains examples that illustrate the 
effects of some of these factors on the observed spectrum. 

The word “spectrum” is used here to mean “Fourier transform.” Because all of the analysis 
presented in this document is in the frequency domain, the Fourier transform has already 
been taken as part of the solution. The figures presented here are plots of the magnitude 
of the spectrum (in dB) versus frequency (in Hz). That is, 20 times the logarithm (base 10) 
of the magnitude of a particular electric field component is plotted on the ordinate versus 
frequency on the abscissa. Therefore, if the spectrum of the x-component of electric field 
is being presented, 20 log [Ej^(t>) E*(v)] will be plotted against frequency in hertz (Hz). 

At very high frequency, the spectra oscillate rapidly, causing a problem in the graphical 
presentation of the results. In particular, as the frequency increases, the resolution (i.e., 
size of the scale) required for resolving these oscillations must also increase. If a fixed scale 
size and a fixed sampling interval are held for the machine calculations, aliasing can occur 
because the results are eventually under-sampled. This happens in most of the figures shown 
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here in which data points at frequencies greater than about 1 MHz are shown, and is especially 
evident in figures 9, 1 1, and 14. In general, the high-frequency results presented here accu- 
rately represent the (upper) envelope of the spectrum, but not the location or number of 
nulls or maxima. Although the algorithms developed for this analysis are capable of including 
the location or number of nulls or maxima, the graphs would have been of impractical size. 

The actual spacing of nulls is discussed in Appendix A. 

Observer's Orientation 

When the observer is on the surface, only a z-component of electric field is measured. How- 
ever, when the observer moves above the surface (for example, onto an elevated platform 
or aboard an airplane), other components of the electric field are observed. This effect is 
illustrated in figures 1 1 through 13, in which a vertical element 1 .5 km long with one end at 
the origin is treated for an observer 100 km from the origin but at angles of 65, 40, and 15 
degrees with respect to the vertical (i.e., axis of the filament). The filament and the observer 
were chosen to be coplanar in all of these cases so that only two components of electric field 
would need to be plotted. The curve labeled E(x) corresponds to a radially (i.e., perpendicular 
to the z-axis) directed field component. For comparison, figure 9 shows the observer located 
on the surface (90 degrees from the vertical). Note that, as the observer approaches the 
axis of the filament, the spectrum becomes flatter (the peak near 10 kHz begins to disappear) 
and there is a general shift toward higher frequency and j^arger amplitude. This effect is 
caused by changes in the effective filament length L(1 - £ • Vp) which can differ significantly 
from the filament’s physical length, L, whenever £ • Vp = 1. 

Figures 14 through 16 show the spectrum of the electric field components radiated by a 
horizontal current element 1 .5 km long with its center 4 km above the origin and parallel to 
the x-axis. The observer is again about 100 km away from the filament and lies in a plane 
that contains the filament and the z-axis; however, the observer is oriented so that the position 
vector of the observer (the line from the origin to the observer) makes angles of 1 5 , 40, and 
65 degrees (figures 14 through 16, respectively) with respect to the vertical. 

Because of symmetry, only two components of electric field exist for this particular orienta- 
tion— the X- and z-components. However, if the observer moves out of the x-y plane so that 
the observer’s position vector and the filament are no longer coplanar, all three rectangular 
field components may exist. Figures 1 7 and 1 8 show examples of such cases for a 1 .5-km 
long horizontal current filament centered 4 km above the origin and parallel to the x-axis. 

The observer is about 100 km away. In figure 17, the observer is on the y-z plane at (y = 

90.63 km; z = 42.26 km) so that the position vector of the observer makes an angle of about 
75 degrees with the vertical. In figure 18, the observer is located at (x = 50 km; y = 30 km; 
z = 2 km). 

Effect of Distance 

The effect of the distance of the observer from the filament is to generally decrease the ampli- 
tude of the spectrum measured by the observer. This is illustrated in figures 19 through 21 
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using the exact solution, and in figures 22 through 24 using the approximate solution. The 
calculations were made for a 1 .5 -km long vertical element, with one end at the origin and 
the observer along a line in the z-x plane which makes an angle of about 1.15 degrees with 
the x-axis. Note the decrease in amplitude with an increase in distance and the general 
improvement in the agreement of the exact and approximate solutions as the distance of 
the observer from the filament increases. As this distance increases, the criterion, kp » 1, 
can be satisfied at lower frequencies, extending the range of agreement between the exact 
and approximate solutions to lower frequencies. 

Orientation of the Current Element 

In most of the preceeding examples, the current element has been vertical with one end 
located on the surface. Although this geometry is apropos of nearly vertical return strokes, 
it is not necessarily representative of other elements of the lightning flash, such as the indivi- 
dual steps in the stepped leader or of intercloud strokes. Figures 25 through 27 show spectra 
which are the result of a filament 1 .5 km long with its center 4 km above the surface. The 
observer was located on the surface, 100 km from the origin. The element and observer 
are coplanar, but the angle that the element makes with the vertical is 45 degrees in figure 
25 and -45 degrees in figure 27; in figure 26 the element is vertical. The case of a horizontal 
element has already been given for this geometry in figure 9. Again note the broadening of 
the spectrum^and the shift toward higher frequencies because of changes in the effective 
length, L(1 - C • Vp), or in the orientation of the filament. 

Effect of Velocity of Propagation 

The velocity with which the current wave propagates along the filament has a definite effect 
on the fields radiated by the element. This effect is demonstrated in figures 28 through 30 
in terms of the approximate (i.e., far -field) solution which, unlike the exact solution, applies 
to all velocities of propagation. In each case shown in the figures, the filament is vertical 
with one end at the origin and is 1 .5 km long, and the observer is located on the surface 
10 km from the origin. Figure 28 shows the electric field measured by the observer for the 
case in which the velocity of propagation, v, equals the speed of light in vacuum, c. In figure 
29, V = c/10, and in figure 30, v = c/100. 

As can be seen, the effect of decreasing the velocity of propagation is to shift the peak in 
the spectrum toward lower frequency with only minor changes in the overaU shape. How- 
ever, if becomes so small that the peak is shifted well beyond the knee of f(v) (figure 8), 
a change in the shape of the spectrum can also be expected. 

The magnitude of the spectrum of the electric field as predicted by the approximate solution 
for an observer on the surface is (equation 28a): 

/~ ~ \]u/e ~ 1 - (t • Vp^)2 ^ 

yj • E* (7^,u) sin — (t?-£*Vp^) 

Pc V - & ‘ Vp^) |_ 2 
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In the case of fi • Vp « 1 , which is reasonable for a distant observer and a vertical element, 
the frequency dependence is roughly: I f(v) I sin (tj kL). That is, the spectrum depends on 
a factor that scales in frequency as kq for fixed L and is multiplied by a factor dependent on 
the waveform (figures 28 through 30). 

TRANSFER FUNCTION FOR THE RADIATION PROCESS 

Examination of the solutions for the radiated fields (i.e., equations 1 1 and 15 for the exact 
solution or equations 3 1 for the approximate solution) indicates that they may be written 
in the form: 


E(r^,i') = f(v)h, (Y^,u) 

(33a) 

H(r„,i^) = l(u)\ ir^,u) 

(33b) 


where f(t') is the Fourier transform of the shape of the disturbance that propagates down the 
filament, and h^ is a vector function that depends on the observer-filament geometry 

and on the physical properties of the filament (such as length, orientation, and velocity of 
propagation) but not on the current waveform. Consequently, the radiation problem is con- 
sidered to be a linear system in which f(iO is the input current waveform, (f^^ , v) is the 
transfer function for the “system,” and E(f^ , u) or H(f^ , u) is the output. T^e requirement, 
inherent in the derivations, that the waveform propagate along the channel, is a characteristic 
that belongs to the system (i.e., all f(i’) propagate). 

Examples of the modulus of the vector components of the transfer function are shown in 
figures 3 1 through 35 for a vertical element 1 .5 km long with one end at the origin. In figure 
32, the observer is on the surface, 100 km form the origin. In this case, only the z-component 
of the transfer function is nonzero, and this component is plotted in the figure for both the 
approximate and the exact solutions. In figures 32 through 34, the observer is also 100 km 
from the filament, but it is located so that the position vector of the observer makes angles 
of 65, 40, and 15 degrees, respectively, with respect to the vertical. 

Note that, for high enough frequency, the transfer function becomes oscillatory with a con- 
stant maximum value for both the approximate and the exact solutions; however, the lower 
bound (not shown in the figures) is zero for the approximate solution and is finite for the 
exact solution (Appendix D). It follows that, at sufficiently high frequency, the envelope 
of the spectrum is determined by the current waveform, independent of geometry. At lower 
frequencies, both geometry and the current waveform affect the shape of the spectrum. 

Goddard Space Flight Center 

National Aeronautics and Space Administration 
Greenbelt, Maryland September, 1976 
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Figure 2. Low-frequency behavior of exact and approximate solutions. 
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Figure 3. Low-frequency behavior of exact and approximate solutions 
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Figure 4. Low-frequency behavior of exact and approximate solutions. 
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Figure 5. Low-frequency behavior of exact and approximate solutions 
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Figure 8. Spectrum of current waveform. 



40 1 . I 1 . . . . , 

lOt 1 lot 2 10’ 3 10»4 lots lots 10*7 lOt 8 10l 9 


FREQUENCY (Hz) 

Figure 9. Exact and approximate solutions: three-exponential waveform. 
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Figure 10. Exact and approximate solutions: three-exponential waveform. 



Figure 11, Exact solution: three-exponential waveform. 
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Figure 12. Exact solution: three-exponential waveform. 
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Figure 13. Exact solution: three-exponential waveform. 
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Figure 14. Exact solution: three-exponential waveform. 
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Figure 15. Exact solution: three-exponential waveform. 
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Figure 16. Exact solution: three-exponential waveform. 
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Figure 17. Exact solution: three-exponential waveform. 
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Figure 18. Exact solution: three-exponential waveform. 



Figure 19. Exact solution: three-exponential waveform. 
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Figure 20. Exact solution: three-exponential waveform. 
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Figure 21 . Exact solution: three-exponential waveform. 
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Figure 22. Approximate solution: three-exponential waveform. 
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Figure 23. Approximate solution: three-exponential waveform. 
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Figure 24. Approximate solution: three-exponential waveform. 



lot 1 101 2 lOf 3 lot 4 lOf 5 lot 6 lOt 7 lOt 8 101 9 

FREQUENCY (Hz) 


Figure 25. Exact and approximate solutions: three-exponential waveform. 
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Figure 26. Exact and approximate solutions: three-exponential waveform. 
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Figure 27. Exact and approximate solutions: three-exponential waveform. 
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Figure 28. Approximate solutions: three-exponential waveform, V = C. 
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Figure 29. Approximate solution: three-exponential waveform, V = C/10. 
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Figure 30. Approximate solution: three-exponential waveform, V = C/100. 
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Figure 31. Exact and approximate solutions: three-exponential waveform, 

f {v) = constant. 
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Figure 32. Exact solution: three-exponential waveform, f (u) = constant. 
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Figure 33. Exact solution: three-exponential waveform, f (i^) = constant. 
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Figure 34. Exact solution: three-exponential waveform, f = constant 
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Figure 35. Exact and approximate solutions: three-exponential waveform 



REFERENCES 


1. Taylor, W. L., “Atmospherics and Severe Storms,” Chapter 17 in Remote Sensing of 
the Troposphere, V. E. Derr, ed., USGPO, 1972. 

2. Jones, H. L., “A Sferics Method of Tornado Identification and Tracking,” Bull. Am. 
Meteorol. Soc. (32) 1951, pp. 380-385. 

3. Jones, H. L., and P. Hess “Identification of Tornadoes by Observation of Waveform 
Atmospherics,” Proceedings of the IRE, 40, 1952, pp. 1049-1052. 

4. Uman, M., Chapter 4. McGraw-Hill, New York, 1969. 


37 




APPENDIX A 
CURRENT WAVEFORM 



APPENDIX A 


CURRENT WAVEFORM 


Measurement of the current in return strokes reveals that the current pulse has a sharp rise 
time (on the order of microseconds to peak value) followed by a relatively slow decay (on 
the order of 50 to 100 microseconds), and a small residual current (continuing current) that 
can persist for hundreds of microseconds after the main pulse is over. A combination of 
three exponentials is common in the literature on lightning for modeling this current (Refer- 
ences 1, 2, and 3). Consequently, such a model has been adopted for use here. 


The following model is taken from Reference 1 : 

I(t) = lo (e-“t - e-^M + Ij e-Tt (Ai) 

where, for a typical lightning return stroke, the parameters are: 

a = 2.0 X 10‘‘ (A2) 

^ = 2.0 X 10^ (A3) 

7 = 1.0 X 103 (A4) 

Iq = 30 ka (A5) 

Ij = 2.5 ka (A6) 


The first two exponentials in equation A1 represent the main current pulse, and the third 
term, Ij exp (~7 t), represents the continuing current. The current waveform, as predicted 
by this model, is as shown previously in figure 7. Note that the peak occurs about 13 micro- 
seconds after initiation of the pulse. 

The spectrum of l(t), which is denoted f(»>) in the text, was shown previously in figure 8. 
Twenty times the logarithm (to base 10) of the magnitude of the Fourier transform of I(t) 
is plotted (along the ordinate) versus frequency (along the abscissa). Curve A represents 
I f(j>) I, and curves B, C, and D represent the (magnitude of) Fourier transform of the three 
constituent exponentials: exp(-at), I^exp(-|3t), and Ij expt (- 7 t), respectively. 

This current waveform has been used in all of the computations made in the text. It is rea- 
sonably representative of the actual current in return strokes, and consequently, the results 
presented here should be reasonable in a quantitative sense for the fields radiated by the 
return stroke. The currents carried by other components of the lightning flash (stepped 
leader, dart leader, and intercloud processes) are not as well known and probably differ 
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from I(t), at least in magnitude (i.e., !„) if not in shape. Consequently, the results predicted 
here would probably be less quantitatively correct when applied to flash elements other than 
return strokes. On the other hand, the analysis is independent of I(t), and, after sufficient 
information is obtained to model I(t) for a given process, the formalism presented here can 
be adopted without change to produce the radiated fields. 

Note that, in the computations of spectra made in the text, it is assumed that the current 
waveform propagates along the channel (filament). That is, the current is I(t - C • rjv), which 
is a wave that travels in the £ direction with speed, v. 
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APPENDIX B 

COORDINATE TRANSFORMATION 


The objective of this appendix is to present a rotation of coordinates that will permit the 
z-axis of the reference coordinate system (x, y, z) to be aligned with an arbitrary unit vector, 
2. One way to accomplish this transformation is to first rotate in the counterclockwise direc- 
tion about the original z-axis through an angle, d = tan'* [ 2^ /2„ ] , until the x-axis of the 

A ^ y 

original system and the projection of 2 on the x-y plane are colinear. The matrix for this 
rotation is: 



0 

sin 6 

cos 0 

0 

cos 6 

-sin 6 

1 

0 

0 


(Bl) 


Next, a counterclockwise rotation is made through the angle, 


tan' 


V X y 


about the new y-axis until the z-axis is aligned with (colinear) the unit vector, The matrix 
for this rotation is: 



- sin 0 0 cos 0 

0 1 0 

cos 0 0 sin 0 


(B2) 


45 





The total rotation matrix is [ R ] = [ R . ] [ ] . Thus, denoting the new coordinates by 


- x"‘ 


"-sin 0 

cos 0sin 6 

cos 0 cos 0" 


" X " 


y" 

== 

0 

cos 6 

-sin0 


y 

(B3) 

. z". 


^COS 0 

sin 0 sin 0 

sin 0 cos 0 , 


_ z ^ 



The inverse matrix is the transpose of [ R ] . 

Of special interest is the transformation of the z-coordinate, which takes the form, 
z “ fi • r, where r = (xx + yy + zz). 
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COMPLETE SOLUTION IN THE FRAUNHOFER REGION 


When p. »L, p. may be expanded in a power series about the distance from the observer to 
the center of the current filament. Keeping only the first few terms, 


Pi 


'K -^c) + -^)l 

I + I — r] I cos + — 

O C Cl 2 



(Cl) 


where 


cos 





and by noting that 


r 


C 


La t + for i 
2 I — for i 


a 

b 


(C2) 


(C3) 


the following form is obtained; 

cos = ±$ • Vp (C4) 

c c 

where p^, = IT^ -7^ I is the distance from the observer to the center of the filament. 

Now, assume that is sufficiently greater that L so that P^t^P^ in all terms in the expres- 
sion for electric field except perhaps where it is multiplied by k, as in the exponentials. 
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expO'k p^); and assume that kL^ / p^ is sufficiently small that the last (and all subsequent) 
terms in the expansion in equation Cl are negligible even in the exponentials. Then, 


e^ 

4 7TPj 


e^ 

4 7tp^ 


exp [ ± jk 


L 

2 



(C5) 


where the (+) pertains for i = a and the (-) for i = b. Finally, assume that k p^ is sufficiently 
large that only terms of lowest order in 1/k p^ need to be considered. Then, only the radia- 
tion terms in e (Pj) are significant and: 


eJk Pi _ 

4 7T Pj ‘ 4 7T p^ 


kL A 

(p^)exp [ ±j — (1 - £ • Vp^) ] (C6) 


where the (+) is associated with i = b, and the (-) is associated with i = a, and 


A 


A 


(P.) = <!)“• r. 


1 + £ • Vp^ 
1 - ^ • Vp^ 


^(pj 


(C7) 




The following relationship for £ • r. has been employed: 


^ • rj = $ • ( rj - r^ + r^) 


A _ L t— for i = a 
£ “ r^ ± — < 

2 1+ for i = b 


(C8) 


A similar result holds for the image terms. By combining the results into an expression for 
E(r^ , v), the following form is obtained for equation 1 1 : 


E ( r^ , yflx/e f (i^) 


eJkPc A 

e„ (pjsin [ i/ 2 kL(l - £ 

4 7TP = 


VPc) ] - 


eJ^^Pc 
4 7tp; 


e'^ (ppsin [ i/ 2 kL(l 


Vp;) 


(C9) 
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where p' is the distance from the observer to the center of the image filament. Comparison 
of equation C9 with the solution obtained by making the far-field approximations before 
integrating (equations 31a and 31b) shows that they are identical. That is, the exact solution 
and the approximate solution converge to the same result in the fraunhofer region. 
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FREQUENCY DEPENDENCE 


This appendix examines the frequency dependence of both the complete and far-field solu- 
tions in the special case of an observer on the surface. Because the location of the observer 
would not be expected to qualitatively effect observations made about frequency dependence 
and because the solutions are easier to handle than otherwise when the observer is located 
on the surface, this case has been chosen for analysis. 

Consider first the complete solution, and examine the magnitude of the spectrum of the 
electric field. Then, as shown in equation E8, 


yjE(r^,v) • E* (r^, ='^n/e I f (»») I | la P + I/? P — 2Re(aj3*) 


cos k (L + Pj, - Pg) - 2 Im (a j3*) sin k (L + p^ ~ } 


’/2 


(Dl) 


where 


a = 



A, . 

iTT 


(D2) 


A 


/3 = 


1 + fi • Vp. 


2 7tp^ 


1 - 




A 


A . . VPu 


e (Ph ) * z + J 


kPh 


(D3) 


and where the relationship, £ • (r^ - r^) = L, has been employed. When the frequency is 
large enough for the imaginary components of a and ^ to be negligible. 


Vi(r- 


E ( r^ , v) • E * ( r^ , i^) ~ 


(D4) 


■^p/e I f (i>) I ■ a^ + — 2 a ^ cos k (L + p^, — p^) 
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This radical oscillates between the extremes of I a + jS I and I a -/3 I , and, consequently, 
the mean behavior at large frequencies is determined by I fCi') 1 , the magnitude of the spec- 
trum of the current waveform. On the other hand, at low frequencies, the radical in equa- 
tion D1 is singular because of the 1/kp dependence of the nonradiation terms, Im a and 
Im ^ . This fact is easily seen by rewriting the radical in the following form: 


V ~ I [ (R® + ~ 2 (Re a) (Re |3) cos k (L + — Pg) ] + 

+ [ (Im a)2 + (ImjS)^ - 2 (Im a) (Im j3) cos k (L + P^ - P^) ] (D5) 

— 2 Im (a j3*) sin k (L + P|, — P^) | 


In the limit of small k, equation D5 becomes: 

~ -|[Reo: — Rej3]^ + [Ima — Imj3]^ — 

- 2 Im (a )3*) sin k (L + p^ ~ 


(D6) 


The first and last terms in equation D6 are independent of k, the limit of small k, and the 
remaining factor, [ Im a - Im /J ] ^ , depends on (1 /k p)^ . Hence, as k p becomes small, the 
radical becomes large, and the spectrum varies as l f(v) I I Im a - Im )3 I . On the other hand, 
if only the radiation terms are considered, the radical approaches the constant value, I Re a 
-Re ^ I , for small k. Consequently, the spectrum that is attributable to only the radiation 
terms varies as I f(p) I in both the large and small frequency extremes. At small frequencies, 
the spectn^ attributable to only the radiation components approaches a limit that is propor- 
tional to I f(0) I , but the spectrum attributable to the complete solution blows up as 1/k 
because of the singularity in the nonradiation fields at zero frequency. 

Next, consider the frequency dependence of the far-field solutions. As shown in equation 
Ell for the far-field solution observed at the surface, 


Vf 


E(T^,p) 




(D7) 


yftx/e I f (i^) 1 


le(p,) • ^1 

2wp^ 




kL sine [ Vi kL (1 
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For large frequencies, the right-hand side of equation D7 oscillates from a maximum of 


2 y[ix/e 


I W I 


l^(Pc) • Z I 

2irp^ 


[ 1 - ^ ♦ VpJ 


to a minimum of zero. To a first approximation, this is the same as the high-frequency be- 
havior of the complete solution shown in Appendix C. However, at low frequencies, this 
solution depends on k I f(i^) I , in contrast to the complete solution which appears as 
I f(i^) I /k if the nonradiation terms are kept or as IF (i^) I if only the radiation terms are 
considered. Examples of the spectrum predicted by the far-field solution with the exponen- 
tial current waveform (Appendix F) were plotted previously in figures 2 through 6. Both 
the radiation terms of the complete solution and the far-field solution are plotted together 
for the exponential current waveform. Note the drop in the spectrum as 1/k in the far-field 
solution at low frequencies and note that both the complete and far-field solutions approxi- 
mate each other at high frequency. 
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APPENDIX E 

FIELDS ON THE SURFACE 


When the observer is on the surface (i.e,, = 0), a symmetry exists at the observation point 

with respect to the source and image that results in simplification of the solutions. In par- 
ticular, when = 0, the following relationships are obtained for i = a, b: 


Pi = p; 

(El) 


^ • Vpj = • Vp.' 

(E2) 


e • r, = £' . r; 

(E3) 

and, because 

Vpj' = VPj - 2 ( Vp. • z) z 

(E4) 


= $- 2($ • $)z 

(E5) 

then 

^-t' = 2(t-z)z 

(E6) 


Vpj - Vpj' = 2(Vp. • $)$ 

(E7) 


By substituting these relationships into equations 1 1 and 15 of the text, the following forms 
are obtained for the complete solution as measured by an observer located on the surface: 


^ ^ 

E(t v)=z \ ix/e T(p^) • z - 

( 2 7TPb 


eJkPa 


2 7T 


e(p J • z > (E8) 


= fiv) 


i eJkPb _ eJkPa _ ) 

^ 2 n ° •’ 2ir ^ ^ 


(E9) 
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where 


ho (Pi) = 

A _ I 1 + t • V/5. [ (Vp. • X ^ - (1^ • ^) V/). X ^ 

ejkK • rj ^ — 1^=^ - - 

"V 1 - 8 ' • Vp| V 1 - (^ * Vpj)2 

Also, substituting the foregoing relationships into equations 31a and 31b of the text, the 
following form is obtained for the far-field solution as seen by an observer located on the 
surface: 


(ElO) 

] 


E(T„,v) = 

^ I ~ eJhPc 

2j z f (^) 7 

4irp^ 

~ ~ eJhPc 

H(7^,«^) = 2jf0^) 

O 4 jrp^ 

Note that, in both equations El 1 and El 2, the electric field intensity is z-directed and that 
h^(p.) • z = 0 for i = a, b, c. Thus, the electric field is perpendicular to the surface, and the 
magnetic field is tangent to the surface as is required at a perfectly conducting boundary. 


(Ell) 

(p^) • ^sin(>/ 2 kL(l - • Vp^)] 

Eq (Pc> sin [ HkLCl • Vp^) ] (E12) 
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DISTRIBUTION OF MINIMA FOR THE 
APPROXIMATE SOLUTION 


The approximate solution for the P* component of the E field can be written as 


E, - (e; - Ef 

( 1 • li - (i ■ 

= K < 

f (t? - K • Vp^ p^ 

1 • [$'-($' • Vp')Vp'] ^ • >■;)) 

c c 

(7j - ^ • vpp p; 


(Fl) 

kL , A 

sin — - (t 7 — C • S/p ) 

2 ® 

sin (t? — fi • V Pg) 


Zeros of E. occur when E.^ = , the most important special case of which is 

E.^ = Ej* = 0, or equivalently, 

>/ 2 kL(t? - • Vp^) = nir 

n, m non zero integers (F2) 

Vi kL (tj — 8 ' • Vp|,) = m ir 

Equation F2 states that a zero of E. occurs when the arguments of the sine terms in equa- 
tion Fl are equal to integer multiples of tt. Alternatively, these conditions may be written 
as 


( 14kL(T? — ^ P^) = nir 

( i/akL(^' • Vp; - ^ • Vp^) = Ktt 


n, £ integers ni=Q 


(F3) 


Using V = ck/27r, the frequencies for which E. = 0 are given by 

c £ 


cn 


V = 


L(t? - COS0) L(COS0’ - cos^) 
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where 


_ A _ 

COS 0 = X • Vp^ 

cos 0' = $ ' • Vp^ 


The geometrical constraint, implicit in the foregoing equations, is that, 

£ _ (cos 0 cos 0 ) 
n (t? — cos 0) 

To consider the simplest case fir^t, let the observer be located on the perfectly conducting 
plane, (Z = 0). Then, £ • Vp^ = £' • Vp' , and conditions of equation F2 are identical 
(i.e., m = n). The frequency at which first goes to zero is 

_ c 

^ L (t? — cos 0) 

Because the zeros are equally spaced, the zero density for a fixed geometry is constant and 
is given by 

d = L(t? — cos 0)/c 

In summary, as L ( t? - cos <p) increases, the first zero shifts toward lower frequencies and 
the zero density increases and vice versa. Finally, note that, as 7?^ 1 {y c), the zero density 
and first zero become increasingly sensitive to changes in the angle. 

For an arbitrary placement of the observer, the analysis is more involved. However, a few 
qualitative statements can be made: 

(a) For an observer above the conducting plane, the probability of satisfying the 
geometrical constraints imposed by equation F2 is zero. Consequently, the probability of 
E. = 0, under conditions (i), is zero. 

(b) To investigate the minima of E. , it is helpful to make use of the following 
definitions: 

A _ 

o)j = OJ (Pc + 1? (£ • r^)) 

COq = (^1 + 

CO 2 = ^ (P'c + vd' • r[,)) 

=> CO2 — Wq “ ^0 ~ ~ Aco/C 
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Then, I Ej I can be expressed in the form: 


kL 


kL 


I E. I = I cos A k [Cj sin (t? — cos 0) — Cf sin (t] — cos 0 ) ] — 


kL 


kL 


— sin Ak [ Cj sin (tj — cos </>) + Cj sin (rj — cos 0') ] I 


(c) Special cases: 

For Cj/C| » 1 , I E. I is minimized at frequencies b V 2 kL (17 — cos 0) *=» 2 ir n. 

For C./C| « 1 , I E. I is minimized at frequencies b ¥2 kL (77 — cos 0') ** 2 tt m. 

For observation points very close to the conducting plane, Ak«l , C./C| « 1 , strong 
minima are expected at the frequencies 

= cn/L (t 7 — cos 0) 

(d) When none of the foregoing conditions hold, the following procedure seems to 
give reliable results: 

= cm/L (t 7 — cos 0') 

= cn/L (t 7 — cos 0) 


An estimate of the first frequency at which E. has a well-defined minima is then given by 

- 

'"‘=1 Cj + c: Cj + Cj 

After V is computed, the density of minima and all subsequent minima can be found. 

For example, consider the spectrum of the field for the case of a 50-m vertical element as 
seen by an observer situated at lO'* ( V 0. ^ (shown previously in figure 35). 

Then, 

V ~ 20.8 MHz 

V ~ 41.6 MHz 

V ~ 62.4 MHz 
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~ 83.2 MHz 


y ~ 104 MHz 

<=5 


These values compare favorably with the values obtained by a high-resolution plot of equa- 
tion FI for this case; i.e., 


(21, 42, 61, 81, 102) MHz 


(e) From the preceding discussion, it is evident that the relative magnitudes of 
C, C' affect the location of the minima. In general. 


r - c 


q + q 


(i 4 j) 


Therefore, a possibility exists that the minima of the various field components will occur at 
significantly different frequencies. Examples of this effect were shown previously in figures 
17 and 18. 
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A(r^,u) = fxf(,u) j e 
filament 


/ [^] 


A A 

£ - £' dr' 


4 7T R 4 7T R 


(18) 


E ( = jk yiti/e f (i^) J 

filament 


r [£•'?'! 

J eJ2w L-7-_ 


$ + — • V ) V 


gjkR 

47tR 


^ ' + — • V ) V 

k^ 


ejkRj 

47t R, 


dr 


(1'9) 


R = lr„ - F' I = l(f^ - r^) + (r, - r') 


- rj2 + ir^ _ r ' [2 + 2( _ r^) . ( - r ') 


= I r — r I + 1 r„ — r ' I cos <i> + — 

o c c 1 


lr„ - f' |2‘ 


IF„ - F, 

o c 


( 20 ) 


cos 0 = 


( ’'o - ^c) • ( r, - r ') 


I r — I 1 — r 

o c c 


( 21 ) 


kL^/Pe « 1 
kL » 1 


(22a) 

(22b) 


E(r^,i^) = jk-^Ji/e 7{y) 


^ + — (^ • V )V 
k^ 


4 7Tp^ 


Hf) - 


(^' • V)V 


-2 — ^ I K 

4 7TP^ 


(23) 
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